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' Abstract. In this paper the N—S asymmetry of the soft X-ray flare index (_F"/sxfl) 

, during the solar cycles 21, 22 and 23 has been analyzed. The results show the 

existence of a real N — S asymmetry which is strengthened during solar minimum. 
The slope of the regression lines fitted to the daily values of asymmetry time series 
. has been found to be negative in all the three cycles. The yearly asymmetry curve 

can be fitted by a sinusoidal function with a period of eleven years. The power 
spectral analysis of daily asymmetry time series reveals the significant periods of 
around 28.26 days, 550.73 days and 3.72 years. 
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■ 1. Introduction 

a^ 

' It is well known that many types of solar phenomena exhibit some 

. form of north-south (N— S) asymmetry. Bell (1962) found a long term 

asymmetry in the sunspot area data for cycle 8 through 18. Roy (1977) 
' studied the N— S distribution of major solar flares from 1955 to 1974, 

covering cycles 19 and 20, and found an asymmetry in favour of north- 



O I ern hemisphere which is more pronounced during the minimum of solar 

' cycle. Knoska (1985) studied the asymmetry of flare activity during 

, the years 1937—1978 using the Ha flare index {FIho) introduced by 

Kleczek (1952, 1953) and found no unique relationship between the 
asymmetry of flare activity and eleven year solar cycle. Vizoso and 
\ Ballester (1987) have studied the sudden disappearance (SD) of solar 

■ prominence during cycles 18—21 and found that the N— S asymmetry 

curve is not in phase with the solar cycle and peaks at or around the 
solar minimum and that the asymmetry changes its sign during the 
maximum, at the time of reversal of the Sun's general magnetic field. 
Verma (1993) reported a study of various solar phenomena occurring 
in both northern and southern hemisphere of the Sun for solar cy- 
cles 8—22. Joshi (1995); Li, Schimider, and Li (1998) and Atag and 
Ozgiic (1998) studied the N— S asymmetry during the solar cycle 22 
with different manifestations of solar activity. In all these studies the 
asymmetry during this cycle was found to be in favour of southern 
hemisphere. Atag and Ozgiig (2001) studied the N— S asymmetry of 
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FiHa during the rising phase of solar cycle 23 (1996—1999) and found 
that the activity prefers the northern hemisphere in general during the 
present cycle. Recently Li et al. (2002) have presented the dominant 
hemisphere of solar activity in each of solar cycle 12 to 22 by calculating 
the actual probability of the hemisphere distribution of sunspot groups, 
sunspot areas, the relative number of sunspots, flare index and number 
of SDs. 

Studies of periodicity in time series of N— S asymmetry data have 
also been made by several authors. Vizoso and Ballester (1989) found 
a significant peak of around 3.1—3.2 years in the power spectrum per- 
formed with the values of flare number and flare index N— S asymmetry. 
Carbonell, Oliver, and Ballester (1993) and Oliver and Ballester (1994) 
studied the periodicity in the N— S asymmetry of sunspot area from 
1874—1989 and 1874—1993 respectively. They obtained statistically 
significant period of 163.93 rotations (12.1 years). Ozgiig and Atag 
(1996) performed the power spectral analysis of N— S asymmetry time 
series of flare index values of solar cycle 22 and confirmed the 25.5 day 
fundamental period of the sun which was discovered by Bai and Sturrok 
(1991). Taking all this into account, here we have made an attempt 
to investigate the N— S asymmetry of soft X-ray flare index (FIsxr), 
which is based on the continuous record of SXR flares observed by 
GOES during solar cycles 21, 22 and 23. The daily values of FIsxR 
include the contribution of SXR flares of different classes (B to X) and 
are suitable for a short and intermediate term solar activity analysis. 
In the paper we have also searched for the periodicity in the daily 
asymmetry time series of FIsxR for the period 1975—2003. 



2. The Soft X-ray Flare Index 

For the present analysis, first we have calculated daily FIsxR by using 
the data of SXR flares observed by GOES in 1—8 A wavelength band. 
The data for the time span of 01 September 1975 to 30 June 2003 has 
been downloaded from NGDC's anonymous ftp server: ftp: / /ftp. ngdc.no] 
aa.gov/STP/ SOLAR_DATA/SOLAR_FLARES/XRAY_FLARES. Dur- 
ing this period the occurrence of 56,072 SXR flares are reported. 

The FlgxR is introduced by Antalova (1996) by weighing the SXR 
flares of classes C, M and X as 1, 10 and 100 respectively (in units 
of 10~^ Wm~'^). In this manner, for example, the SXR flare index 
for individual flares of class C7.3, M2.9 and X4.8 is 7.3, 29.0 and 480 
respectively (See also Landi et al., 1998). Here, we have also included 
the contribution of flares of class B in the calculation of daily FIsxR- 
Suppose niB, rnc, tum and mx are the digit multipliers for the flares 
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Figure 1. Plot of the monthly values of FIsxr for the period 1975—2003. Dashed 
line shows the smooth spline curve. 



of class B, C, M and X respectively, then the expression to calculate 
daily FIsxr can be written in the following way 



FIsxR = 0.1 X ^ me + 1.0 x ^ mc + 10.0 x ^ ttim + 100.0 x ^ mx 

i j k I 

(1) 

where i, j, k and / are the number of flares of class B, C, M and X per 
day respectively. 

Figure 1 shows the variation of monthly values of FIsxR for the pe- 
riod of 1975—2003 covering solar cycles 21, 22 and the rising and maxi- 
mum phase of cycle 23. We have compared FIsxr with FIho- For this 
a correlation plot of these activity indicators is given in Figure 2. The 
FIho values are available for general use in Kandilli observatory's and 
NGDC's anonymous ftp servers: ftp://ftp.koeri.boun.edu.tr/pub/astr- 
onomy/flareJndex and ft p://ftp ngdc.noaa.gov/STP/S OLA R_DATA7 
SOLAR_FLARES /INDEX^ 
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Figure 2. Scatter diagram showing the correlation between the monthly values of 
FIsxR and Flua- Here r indicates the correlation coefficient. 



3. The Analysis of Asymmetry Time Series 

The N— S asymmetry for the FIsxR is defined as 



AsXR = 



FIn - FIs 
FIn + FIs' 



(2) 



where FIn and FIs stands for the daily FIsxR in the northern and 
southern solar hemisphere respectively. Thus, if AsxR >0, the activity 
in the northern hemisphere dominates and if AgxR <0, the reverse is 
true. In the GOES flare list the heliographic latitudes of all the flares 
are not known and SXR flares of class > B 1.0 are not occurred daily. 
Therefore such data were excluded from our study. Due to these reasons 
the above expression gives us an asymmetry time series composed of 
6952 values. In Figure 3 the plot of daily asymmetry time series is 
presented. To show the reality of the variations of asymmetry time 
scries we have fitted a straight line to the daily values of AsxR for each 
cycle separately. 

To investigate to what extent the asymmetry is real wc have followed 
the method of Letfus (1960) in which we can define the asymmetry of 
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Figure 3. The plot of the daily values of asymmetry time series during solar cycles 
21, 22 and 23. The straight line corresponds to the fit of first order polynomial. 

random distribution on the solar disk as 

l^AsxR = ± ^===1=== , (3) 

which depends upon the vahies of FlgxR in the northern and southern 
hemisphere respectively. To verify the reliability of calculated asymme- 
try values, test is applied with 

^ 2{FIn - FIs) ^ V2AsxR 
^ ^2{FIn + FIs) ^AsxR ■ ^ ' 

Thus for AsxR < AAsxR, ^AsxR < AsxR < 2AAsxR and AsxR > 
2AAsxR, the probability that N— S asymmetry exceeds the dispersion 
value is p < 84 %, 84 % < p < 99.5 % and p > 99.5 % respectively. Here 
the first, second and third limits imply for the statistically insignificant, 
significant and highly significant values respectively. Using the method 
of Letfus we have found that out of total 6952 cases, p > 99.5 % in 4789 
cases; 84 % < p < 99.5 % in 1237 cases and p < 84 % in 926 cases, which 
means that, in general the N— S asymmetry is a real phenomenon and 
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Figure 4- Yearly N — S asymmetry index of FIsxr during 1976—2003. Highly signif- 
icant asymmetry values with p> 99.5 % are marked with circles, otherwise squares 
are drawn. Dashed line corresponds to the fitted sinusoidal curve with a period of 
eleven years. Downward arrows indicate solar activity minima. 

not due to random fluctuations. We have also studied the asymmetry 
of FlgxR year wise (see Figure 4) . In this case we have found that out 
of 28 cases, the data is highly significant (p > 99.5 %) in 27 cases. 

We are also interested to see any possible connection of the asymmet- 
ric behaviour of FIsxr with a period of eleven years. For this, in Figure 
4, a comparison between the original asymmetry curve and a fitted 
sinusoidal function is also presented. The fitted sinusoidal function is 
of the form 

y{x) = C + Asm2Trf{t-tQ), (5) 

where C is a constant term and A is the amplitude, to has been assigned 
a value 1976. The frequency / has been given a fixed value correspond- 
ing to eleven year period. Thus / = 1/11 or 0.091 years The values 
of C and A come out to be —0.22 and 0.35 respectively. We find that 
there is a good agreement between original curve of asymmetry index 
and the fitted sinusoidal curve. 
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4. Power Spectra of Asymmetry Time Series 

Using the Lomb-Scargle periodogram method (Lomb, 1976; Scaxgle, 
1982; Horne and Bahunas, 1986) for the daily asymmetry time series, 
which is plotted in Figure 3, we have performed power spectra analysis. 
Figures 5a, 6a and 7a show the normalized power spectra of asymme- 
try time series during 1975—2003 for the three frequency ranges of 
116.0—463.0, 16.5—116.0 and 7.7—16.5 nHz which corresponds to the 
period intervals of 100-25, 700-100 and 1500-700 days respectively. 
The FIsxR is not independent. Therefore the probability P of the 
power density at a given frequency being greater than K by chance is 
given by 

P{z >K) = exp{-K/k), (6) 

where the normalization factor k, which is due to event correlation, can 
be determined empirically (Bai and Oliver, 1990). 

Figures 5b, 6b and 7b show the distribution of the power values 
corresponding to the normalized power spectra shown in Figures 5a, 
6a and 7a respectively. The vertical axis shows the cumulative number 
of frequencies for which the power exceeds a certain value. For example, 
in the power spectrum for 116.0—463.0 nHz range shown in Figure 5a 
we have 1222 frequencies. For all these frequencies the power exceeds 
zero; thus, we have a point at (X = 0, Y = 1222). At only one frequency 
(409.5 nHz, which is equivalent to 28.26 days) the power was 30.02, its 
maximum value. For lower values of power, the distribution can be well 
fitted by the equation y = 1226 exp(-x/3.7), as expected from equation 
(6). Thus, we normalize the power spectrum by dividing the powers by 
3.7 to obtain Figure 5a. For other cases we use the same procedure for 
normalization. 

Once the power spectrum has been normalized properly, we can use 
the 'False alarm probability' (FAP) formula to estimate the statisti- 
cal significance of a peak in the power spectrum. It is given by the 
expression 

F = 1 - [1 - exp{-Zm)f, (7) 

where is the height of the peak in the normalized power spectrum 
and N is the number of independent frequencies. Therefore, for exam- 
ple, in the case of power spectrum shown in Figure 5a, if we substitute 
Zjn=8.l3 and N=1222 (since we searched 1222 frequencies with 0.284 
nHz intervals) in equation (7) we get F=0.30, i.e., the probability to 
obtain such a high peak at 409.5 nHz (28.26 days) by chance is about 
30%. The same analysis has been applied to the power spectra in other 
two frequency ranges. In the frequency interval of 16.5—116.0 nHz and 
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Figure 5. (a) Normalized power spectrum of the daily asymmetry time series for the 
frequency interval of 116.0—463.0 nHz (100—25 days), (b) Scargle power distribution 
corresponding to Figure 5a. The vertical axis is the number of frequencies for which 
power exceeds x. The straight line is the fit to the points for lower values of power. 



paper.tex; 5/02/2008; 13:23; p. 8 



9 



Period (Days) 



578,7 289.4 193.0 144.7 115.7 




20.0 40.0 60.0 80.0 100.0 

Frequency (nHz) 




10 20 30 40 

Power 



Figure 6. Same as Figure 5, but for the frequency interval of 16.5—116.0 nHz 
(700-100 days). 



7.7—16.5 nHz we get important peaks at 21.02 nHz (550.73 days) with 
a FAP =7% and 8.5 nHz (3.72 years) with a FAP=20.4% respectively 
To confirm the peaks of 28.26 days, 550.73 days and 3.72 years 
in Figures 5a, 6a and 7a respectively are not due to aliasing, we have 
removed sine curves of these periods from the original time scries. Then, 
once again, we have computed the periodogram (see Figure 8). 
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Figure 7. Same as Figure 5, but for the frequency interval of 7.7—16.5 nHz 
(1500-700 days). 



5. Discussions and Conclusions 



In this paper we have calculated the daily FIsxR for solar cycles 21 to 
23. The plot (Figure 1) of monthly values of FIsxR indicates a lower 
level of activity during solar cycle 23 compared to previous two cycles. 
This has been pointed out earlier by Ozgiic et al. (2002) in his study 
during the rising phase of solar cycle 23. The scatter diagram (Figure 



paper.tex; 5/02/2008; 13:23; p. 10 




11 



a 



135,0 180,0 225,0 270,0 315,0 360,0 405,0 450,0 



^ 8 





20.0 



40,0 



60,0 



80,0 100,0 



Frequency (nHz 



c 




Figure 8. Pcriodogram of asymmetry time series, for the frequency intervals of (a) 
116.0-463.0 nHz (b) 16.5-116.0 nHz and (c) 7.7-16.5 nHz, after subtraction of a 
sine curve with frequency 409.5 nHz (28.26 days), 21.02 nHz (550.73 days) and 8.5 
nHz (3.72 years) respectively. 
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2) between FIsxR and Flua shows that there exists a good correlation 
between these two solar activity indices. Atac and Ozgiic (1989, 1998, 
2003) in their various papers have made extensive statistical analysis of 
Flffa, finding its good correlation with a number of solar indices that 
arise under different physical conditions. Therefore, a good correlation 
of FIsxR with FIho also makes the former an interesting parameter 
which is good enough to describe the solar activity. 

We have found that there exists a real N— S asymmetry in FIsxR 
which is not due to random fluctuations. The fitted straight line to 
asymmetry time scries (Figure 3) for the cycles 21 and 23 shows that 
the activity in the northern hemisphere is more important during the 
ascending branch of cycle whereas during the descending branch the 
activity becomes more important in the southern hemisphere. In the 
cycle 22 there was a southern dominance throughout. The slope of 
regression lines for all the three cycles is found to be negative. Vizoso 
and Ballester (1990), in his study of N— S asymmetry of sunspot area 
during solar cycles 12—20 and sunspot number for cycle 21, found that 
the slope of regression lines fitted to the yearly values of N— S asym- 
metry change its sign each four cycles suggesting a long term periodic 
behaviour in the N— S asymmetry of around eight cycles. They also 
found that the slope of regression lines has changed from positive to 
negative values after cycle 19. Atag and Ozgiig (1996) made similar 
conclusions with the N— S asymmetry of Flna during cycles 17—22. 
Our study confirms and extends these earlier results. Moreover, Atag 
and Ozgiig (1996) suggested that, if the period of eight cycles is real, 
the behaviour of asymmetry during cycles 20—22 should be maintained 
in cycle 23, which has also been confirmed by the present study. 

On the other hand, if we consider the plot of yearly N— S asymmetry 
index (Figure 4) the change in the predominance of activity in the two 
hemisphere can be studied for every year starting from 1976. Initially 
during the solar cycle 21, the northern hemisphere was more active but 
after 1980 it was shifted slightly towards south. In the year 1983 the 
activity moved towards south strongly and prevailed there during most 
of the solar cycle 22. In the solar cycle 23, we observe a strong peak in 
southern hemisphere during the minimum phase in the year 1996 and 
then the activity is moved towards the north. In the years 1998, 1999 
and 2000 the northern hemisphere was dominated. These results are in 
agreement with the work done by Temmer et al. (2001), who analyzed 
the N— S asymmetry of Ha flares from 1975 to 1999. The preference for 
northern hemisphere during the rising phase of cycle 23 is also reported 
by Atag and Ozgiig (2001) with the data of Flua- Here we also notice 
that, during cycle 23, the behaviour of asymmetry is changed in the 
year 2000 and it has dominated towards southern hemisphere in the 
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successive years. The Figure also reveals that the fitted sinusoidal curve 
with a period of eleven years reflects the behaviour of asymmetry in 
a good manner. This suggests that N— S asymmetry of FIsxR has a 
periodic behaviour shifted in phase with respect to solar cycle. The 
shape of the curve also indicates that the asymmetry has peaked at 
or around the minimum of solar activity. This resuh is complementary 
to the study of N— S asymmetry made by several authors (Swinson, 
Koyama and Saito, 1986; Vizoso and Ballester, 1990; Joshi, 1995; Atag 
and Ozgiig, 1996). 

Finally, wc have studied periodicities in the daily asymmetry time sc- 
ries for three frequency ranges. In the power spectrum for the frequency 
range of 116.0—463.0 nHz (Figure 5), we get a significant peak at 409.5 
nHz (28.26 days). It is close to the 27 days solar rotational period which 
is well known in different solar indices. Antomicci ct al. (1990) investi- 
gated the rotation of photospheric magnetic field during solar cycle 21, 
and obtained a dominant period of 26.9 days for the northern and 28.1 
days for the southern hemisphere. A similar outcome is obtained by 
Temmer et al. (2002) from the power spectral analysis of daily sunspot 
numbers during 1975—2000 in which they derive a rigid rotation with 
27.0 days for the northern hemisphere, while the southern hemisphere 
reveals a dominant period of 28.2 days. Pap, Tobiska, and Bouwer 
(1990) and Bouwer (1992) reported periods between 26—28 days from 
the power spectra of various solar indices. Joshi (1999) reported peaks 
at 26.5 days and 28.3 days in the power spectrum of solar radio flux at 
10.7 cm. Lou et al. (2003) found a prominent peak of around 28 days 
in the daily averaged Ap index for geomagnetic disturbances, which is 
physically identified with solar rotation. The power spectrum for the 
16.5—116.0 nHz range (Figure 6) shows a significant peak at 21.02 nHz 
(550.73 days). Ichimoto et al. (1985) reported that a peak in the power 
spectra appears near 510—540 days when northern and southern hemi- 
sphere flares are analyzed separately. Bai (1987) found an important 
peak at 552 days in the power spectra of major flares belonging to 
northern hemisphere. A significant peak of around 20 rotations (about 
540 days) is found by Akioka et al. (1987) in the power spectra for TA 
(Sum of maximum areas of sunspot groups per solar rotation) and MA 
(mean area of sunspot groups per solar rotations) during solar cycle 21. 
Ozgiig and Atag (1989) reported a peak at 564 days in the FIho during 
cycle 20 and 21 considered together. Oliver, Carbonell, and Ballester 
(1992) found the periodicities around 540 days in the cycles 12, 14 and 
17 in sunspot areas, while during cycles 18 and 19 it is present, with 
a very high significance, in sunspot and Zurich sunspot number. They 
also found a peak at 528 days in the periodogram for sunspot areas 
during cycles 12 to 21. In the power spectra for the frequency range 
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of 7.7—16.5 nHz (Figure 7), we have a significant peak at 8.5 nHz 
(3.72 years). Rao (1973) studied periodicities in several indices of solar 
activity and found that the index P (average area of all spot groups), Tq 
(average spot group life time) and /o (number of all spot groups formed 
all over the Sun within a unit time) presented periodicities around 3.9, 
3.5 and 3.48 years, respectively. Vizoso and Ballester (1989) found a 
periodicity of around 3.1—3.2 years in the values of Flua and flare 
numbers N— S asymmetry. A periodicity of 3.27 years is reported by 
Vizoso and Ballester (1990) in the N— S asymmetry of sunspot areas. 
In all the three power spectra presented by us, the subtraction of sine 
curve indicates that the peaks are probably real and not spurious or due 
to aliasing (Figure 8). We, therefore, see that the periodicities in the 
asymmetry time series of FlgxR is consistent with the earlier results. 



Acknowledgements 

We are thankful to Prof. Ram Sagar for valuable comments and sug- 
gestions. One of the authors (BJ) wishes to thank Mr. Brijesh Kumar 
and Mr. J. C. Pandey for help in data analysis. 



References 

Akioka, M., Kubota, J., Suzuki, M., Ichimoto, K., and Tohmura, I.: 1987, Solar 
Phys. 112, 313. 

Aiitalova, A.: 1996, Contrib. Astron. Ohs. Skalnate Pleso. 26, 98. 

Antonucci, E., Hocksema, J. T., and Schcrrcr, P. H. : 1990, Astrophys. J. 360, 296. 

Ata5, T. and Ozgiig, A.: 1996, Solar Phys. 166, 201. 

Atag, T. and Ozgiig, A.: 1998, Solar Phys. 180, 397. 

Atag, T. and Ozgiig, A.: 2001, Solar Phys. 198, 399. 

Bai, T.: 1987, Astrophys. J. 318, L85. 

Bai, T. and Oliver, E. W..: 1990, Astrophys. J. 363, 299. 

Bai T. and Sturrock, P. A.: 1991, Nature 352, 360. 

Bell, B.: 1962, Smithsonian Contr. Astrophys. 5, 203. 

Bouwer, S. D.: 1992, Solar Phys. 142, 365. 

CarboncU, M., Oliver, R., and Ballester, J. L.: 1993, Astron. Astrophys. 274, 497. 
Horne, J. H. and Baliunas, S. L.: 1986, Astrophys. J. 302, 757. 
Ichimoto, K., Kubota, J., Suzuki, M., Tohmura, I., and Kurokawa, H.: 1985, Nature 
316, 422. 

Joshi, A.: 1995, Solar Phys. 157, 315. 

Joshi, A.: 1999, Solar Phys. 185, 397. 

Klcczek, J.: 1952, Publ. Contr. Astron. No. 22, Prague. 

Kleczek, J.; 1953, Publ. Astrophys. Obs. Czech. Acad. Sci. No. 24, Prague. 

Knoska, S.: 1985, Contrib. Astron. Obs. Skalnate Pleso 13, 217. 



paper.tex; 5/02/2008; 13:23; p. 14 



15 



Landi, R., Moreno, G., Storini, M., and Antalova, A.: 1998, J. Geophys. Res. 103, 

No. A9, 20553. 
Letfus, v.: 1960, Bull. Astron. Inst. Czech. 11, 31. 

Li, K. J., Schmieder, B., and Li, Q. S.: 1998, Astron. Astrophys. 131, 99. 

Li, K. J., Wang, J.X., Xiong, S. Y., Liang, H. F., Yun, H. S., and Gu, X. M.: 2002, 

Astron. Astrophys. 383, 648. 
Lomb, N.: 1976, Astrophys. Space Set. 39, 477. 

Lou, Y., Wang, Y., Fan, Z., Wang, S., and Wang, J.: 2003, Monthly Notices Royal 

Astron. Soc. , In press. 
Oliver, R., Carbonell, M., and Ballester, J.L.: 1992, Solar Phys. 137, 141. 
Oliver, R., Ballester, J.L.: 1994, Solar Phys. 152, 481. 
Ozgiig, A. and Atac, T.: 1989, Solar Phys. 123, 357. 
OzgiiQ, A. and Atag, T.: 1996, Solar Phys. 163, 183. 

Ozgiig, A., Atag, T., and Rybak, J.: 2002, J. Geophys. Res. 107, No. A7, 

10.1029/2001JA009080. 
Ozgiig, A., Atag, T., and Rybak, J.: 2003, Solar Phys. 214, 375. 
Pap, J., Tobiska, W.K., and Bouwer, S. D.: 1990, Solar Phys. 129, 165. 
Rao, R.: 1973, Solar Phys. 29, 47. 
Roy, J. R.: 1977, Solar Phys. 52, 53. 
Scarglc. J. D.: 1982, Astrophys. J. 263. 835. 

Swinson, D. B., Koyama, H, and Saito, T.: 1986, Solar Phys. 106, 35. 

Temmer, M., Vcrouig, A., Hanslmeier, A., Otruba, W., and Messerotti, M..: 2001, 

Astron. Astrophys. 375, 1049. 
Temmer, M., Veronig, A., and Hanslmeier, A..: 2002, Astron. Astrophys. 390, 707. 
Vcrma, V. K.: 1993, Astrophys. J. 403, 797. 
Vizoso, G. and Ballester, J. L.: 1987 Solar Phys. 112, 317. 
Vizoso, G. and Ballester, J. L.: 1989 Solar Phys. 119, 411. 
Vizoso, G. and Ballester, J. L.: 1990 Astron. Astrophys. 229, 540. 

Address for Offprints: Bhuwan Joshi 

State Observatory, 

Manora Peak, Naini Tal-263 129, 

INDIA. 

Phone: +91-05942-235583, 235136 
Fax: -h91-05942-235136 
email: bhuwan@upso.ernet.in 



paper.tex; 5/02/2008; 13:23; p. 15 



paper.tex; 5/02/2008; 13:23; p. 16 



